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Abstract 



We propose to test the homogeneity of a Poisson process observed on a finite interval. 
In this framework, we first provide lower bounds for the uniform separation rates in 
f-H I norm over classical Besov bodies and weak Besov bodies. Surprisingly, the obtained lower 

X/^ • bounds over weak Besov bodies coincide with the minimax estimation rates over such classes. 

Then we construct non asymptotic and nonparametric testing procedures that are adaptive 
in the sense that they achieve, up to a possible logarithmic factor, the optimal uniform 
separation rates over various Besov bodies simultaneously. These procedures are based on 
model selection and thresholding methods. We finally complete our theoretical study with a 
Monte Carlo evaluation of the power of our tests under various alternatives. 
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^ ! 1 Introduction 

a^ 

O ■ Poisson processes have been used for many years to model a great variety of situations: machine 

breakdowns, phone calls... Recently Poisson processes become popular for modeling occurrences 
of words or motifs on the DNA sequence (see Robin, Rodolphe and Schbath [21]) • In this context, 
it is particularly important to be able to detect abnormal behaviors. 
! With such applications in mind, we consider in this paper the question of testing the homogeneity 

of a Poisson process N. Since we can only observe a finite number of points of the process, this 
question has a sense only on a finite interval. For the sake of simplicity, we assume that the 
Poisson process N is observed on the fixed set [0, 1] , and that it has an intensity s with respect 
to some measure /i on [0, 1] with dfi{x) = Ldx. 

Denoting by So the set of constant functions on [0, 1] , our aim is consequently to test the null 
hypothesis {Hq) "s G 5o", against the alternative {Hi) "s 5o". 

This problem of testing the homogeneity of a Poisson process has been widely investigated both 
from a theoretical and practical point of view (see Bain, Engelhardt, and Wright |2j or Cohen and 
Sackrowitz [7] for a survey and Bhattacharjee, Deshpande, and Naik-Nimbalkar [5] for a more 
recent work). In these papers, the alternative intensities are monotonous. Another related topic 
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is the problem of testing the simple hypothesis that a stationary process is a Poisson process 
with a given intensity. We can cite for instance the papers by Fazli and Kutoyants [l2j where the 
alternative is also a Poisson process with a known intensity, Fazli [H] where the alternatives are 
Poisson processes with one-sided parametric intensities, or Dachian and Kutoyants [lO], where 
the alternatives are self-exciting point processes. The paper by Ingster and Kutoyants [H] is the 
closest one to the present work. The alternatives considered by Ingster and Kutoyants are Poisson 
processes with nonparametric intensities in a Sobolev or Besov B2^^{R) ball with 1 < q < +00 
and known smoothness parameter 6. 

However, in some practical cases like the study of occurrences of words or motifs on a DNA 
sequence, such smooth alternatives cannot be considered. The intensity of the Poisson process in 
these cases may burst at a particular position of special interest for the biologist (see Gusto and 
Schbath [H] for more details). The question of testing the homogeneity of a Poisson process then 
becomes "how can we distinguish a Poisson process with constant intensity from a Poisson process 
whose intensity has some small localized spikes?". This question has already been partially 
considered in the seventies in a precursory work by Watson [27]: he proposed a test based on 
the estimation of the Fourier coefficients of the intensity without evaluating the power of the 
resulting procedure. 

In this paper, we focus on constructing adaptive testing procedures i.e. which do not use any 
prior information about the smoothness of the intensity s, but which however have the best 
possible performances (in a minimax sense). 

From a theoretical point of view, we evaluate the performances of the tests in terms of uniform 
separation rates with respect to some prescribed distance d over various classes of functions. 
Given /? g]0, 1[, a class of functions 5i, and a level a test with values in {0,1} (rejecting 
[Hq) when = 1), the uniform separation rate p{(^a-,Si, (3) of over the class Si is defined 
as the smallest positive number p such that the test has an error of second kind at most equal 
to /3 for all alternatives s in Si at an distance p from Sq. More precisely, if denotes the 
distribution of the Poisson process A'^ with intensity s, 



/9($„,5i,/3) = inf<^p>0, sup Ps($a = 0)</?^ (1.1) 

I s£Si,d{s,So)>p I 



= inf<^p>0, inf P^(^>c, = 1) > 1 - /3 ^ . (1.2) 

[ s&Si,d{s,So)>p J 

In view of the practical situations of our interest, we study some classes of alternatives that 
can be very irregular, for instance that can have some localized spikes. We then consider some 
classical Besov bodies and also some spaces that can be viewed as weak versions of these classical 
Besov bodies and that are defined precisely in the following. The interested reader may find in 
Rivoirard [23] some illustrations of functions in weak Besov spaces and how the smoothness 
parameters of the functions govern the proportion and amplitude of their spikes. 
As a first step, we evaluate the best possible value of the uniform separation rate over these 
spaces. In other words, we give a lower bound for 

p(5i,a,/3) =inf/5($,,5i,/3), (1.3) 

where the infimum is taken over all level a tests <&q,, and where Si can be either a Besov body 
or a weak Besov body. This quantity introduced by Baraud [3] as the (a, /3)-minimax rate of 
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testing over Si or the minimax separation rate over Si is a stronger version of the (asymptotic) 
minimax rate of testing usually considered. The key reference for the computation of minimax 
rates of testing in various statistical models is the series of papers due to Ingster [16]. Concerning 
the Poisson model, Ingster and Kutoyants pS] give the minimax rate of testing for Sobolev or 
Besov B^^q{R) balls with I < q < +00 and smoothness parameter S > 0. They find that this 

-25 

rate of testing for the Sobolev or Besov norm or semi-norm is of order L4i+i . Let us note that 
we find here lower bounds for the classical Besov bodies similar to Ingster and Kutoyants'ones. 
Furthermore, our lower bounds for the weak Besov bodies are larger than the ones for classical 
Besov bodies. Alternatives in weak Besov bodies are in fact so irregular that it is as difficult to 
detect them as to estimate them. The problem of estimation in weak Besov spaces is solved by 
using thresholding procedures: indeed the weak Besov spaces are closely related to the maxisets 
of those procedures (see Kerkyacharian and Picard [19] in the Gaussian framework and Reynaud- 
Bouret and Rivoirard [22] in a Poisson model). To our knowledge, no previous results of this 
kind exist for weak Besov bodies in testing problems, even in more classical statistical models, 
like the density model. Despite the similarity of both models, our lower bounds over weak Besov 
bodies cannot however be straightly transposed to the density model since our proofs heavily 
rely on the Poissonian independence properties. 

As a second and main step, we construct non asymptotic level a tests which achieve, up to a 
possible logarithmic factor, the minimax separation rates over many Besov bodies and weak Besov 
bodies simultaneously, whereas using no prior information about the smoothness of the intensity 
s. Our idea here is to combine some model selection methods that are effective for alternatives 
in classical Besov bodies and a thresholding type approach, inspired by the thresholding rules 
used for adaptive estimation in weak Besov bodies. Key tools in the proofs of our results are 
exponential inequalities for U-statistics of order 2 due to Houdre and Reynaud-Bouret [T5] . 
Of course, both model selection and thresholding approaches have already been used to con- 
struct adaptive tests in various statistical models. One can cite among others the papers by 
Spokoiny ([25] and [26J) in Gaussian white noise models or by Baraud, Huet and Laurent [3] in 
a Gaussian regression framework. These papers propose adaptive tests which combine methods 
closely related to both model selection and thresholding ones. As for the density framework, 
adaptive tests were proposed by Ingster [17] or Fromont and Laurent [13], using model selection 
type methods and by Butucea and Tribouley [6] using thresholding type methods. 
The present work is organized as follows. In Section [2l we provide lower bounds for the uniform 
separation rates over various Besov bodies. Our testing procedures are defined in Section [31 and 
their uniform separation rates over Besov bodies are established in Section [H We carry out a 
simulation study in Section [5] to illustrate these theoretical results, and the proofs are postponed 
to the last section. 

2 Lower bounds for the minimax separation rates over Besov 
bodies 

We consider the Poisson process N with intensity s with respect to some measure ^ on [0, 1] , 
with dfi{x) = Ldx. In the following, we assume that s belongs to L^([0, 1]), and (., .), ||.|| and d 
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respectively denote the scalar product 

(/,ff> = / f{x)g{x)dx, 
J[o,i] 

the L^— norm 

ll/f = / f\x)dx, 
J[0,1] 

and the associated distance. 

Let us denote the Haar basis of L^([0, 1]) by {(po, 4>[j,k),j G N, /c G {0, . . . , 2-' — 1}} with 

(j)o{x) = l[Q^l]{x), 

and 

</)(,,fe)(x) = 2^'/2^(2^-3,_^)^ (2.1) 
where il){x) = l[o,i/2[(a^) - l[i/2,i[(a;). 

We set ao = (s, 0o) and for every j G N, /c € {0, . . . , 2-' — 1}, a(j,jt) = (s, </'(j,fc))- 
We can now introduce the Besov bodies defined for 5 > 0, > by 

s > 0, s G L2([0, 1]), s = ao4>o + 

jgN A;=0 

2^-1 1 

VjGN,^a2^.,)<i?22-2^n, (2.2) 
fc=o J 

and more generally for p > 1, i? > and 5 > max(0, 1/p — 1/2), 

{2J-1 
s > 0, s G l2([0, 1]), s = ao0o + X] X] ^{hk)(t>{j,k), 
j&i k=0 



Vj G N, < i?P2'^n^+^i} . (2.3) 

k=0 ) 



As in Reynaud-Bouret and Rivoirard [22], we also introduce some weaker versions of the above 
Besov bodies given for 7 > and i?' > by 

!2J-1 
s > 0, s G l2([0, 1]), s = ao^o + Y,Y1 ^U,k)(t>u,k), 
jgN fe=0 

2^-1 \ 
jr'eN k=0 ) 



Fixing some levels of error a and /? in ]0, 1[, and denoting by L°°(i?") the set of functions 
bounded by R" , our purpose in this section is to find sharp lower bounds for p{B2oa{R) H 
W.,{R') n L°°(i?"),a,/3), where p is defined by (fOll . 

Starting from a general idea developed by Ingster [16], we obtain the following result. 
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Theorem 1. Assume that R > 0, R' > 0, and R" > 2, and fix some levels a and (3 in ]0, 1[ such 

that a + f5< 0.59. 

(i) 7/7 > max(25, 1/2), then 



/ L \ 1 + 27 r 

liminf — p{Bi ^{R) nWJR') n]L'^{R"),a, p) > 0. 
(ii) If5> max(7/2,7/(l + 27)), then 



lim inf L^p{Bi^{R) n VF^(fl') n L°°(fl"), a, /3) > 0. 

L— >+oo ~ ' 

(m) //<5 < 7/(1 + 27) and 7 < 1/2, then 

\i-m-uiiL^p{Bi^{R)r\W.i{R')r\U^{R"),a,P) > 0. 

L— >+oo ~ ' 

Comments. 



1. For the whole set of parameters ((5,7) such that 5 > 7/(1 + 27), we prove in Section [4] that 
these lower bounds are actually sharp. 

2. We have in case (n) lower bounds which coincide with the minimax rates of testing obtained 
by Ingster and Kutoyants pLSj when testing that a periodic Poisson process has a given 
intensity in the Besov spaces B^l^{R). We know (see Ingster [17] or Fromont and Laurent 
[T3] for instance) that such rates can be achieved by some multiple testing procedure based 
on model selection type methods. This is the principle of our first procedure described in 
Section 13. 1[ 

3. We notice that the lower bounds obtained in case (i) are equal to the minimax estima- 
tion rates on the maxisets of the thresholding estimation procedure, namely B'^^^'^'^"'^ (R) H 
W^{R') (see Kerkyacharian and Picard [1^, Rivoirard [23], or Reynaud-Bouret and Rivoirard 
[22] for more details). This means that it is as difficult to test as to estimate over such 
classes of functions, phenomenon which is quite unusual. Since the minimax estimation 
rates on these classes are achieved by thresholding rules, it will be natural to construct 
a testing procedure based on thresholding methods: this is the idea that originated our 
second procedure described in Section 13. 2[ 



3 Two tests of homogeneity 

Let us recall that Sq denotes the set of constant functions on [0, 1] and that we assume that s 
belongs to h^{[0,l]). 

In this section, we construct level a tests of the null hypothesis (Hq) "s G 5o", against the 
alternative (Hi) "s 5o", from the observation of the Poisson process or the points {Xi, I = 
1, . . . , Nl} of the Poisson process. 

We introduce two testing procedures that come from two different statistical approaches. The 
first one originates in general model selection methods, while the second one is closer to the 
thresholding type methods. 
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In order to understand the global ideas of these procedures, let us notice that the squared 
L^— distance (P{s,So) between s and the set of constant functions Sq can be rewritten as 



(fis,So) = / \s{x)- / s{y)dy] dx 
J\o.i] \ J\o.i] 



[0,1] \ J [0,1] 

I ii2 2 
|s|| -Uq, 



where oq = (s, (po), and for all A G Aoo = {{j, k), j G N, k £ {0, ... ,2^ - 1}}, ax = (s, (px). 
For every A G Aqo, oa can be estimated by 

orx = \ [ 4>x{x)dN^, 
^ J [0,1] 

which is also equal to 

1 

OA = - ^(pxiXl). 

1=1 

From this variable, we deduce an unbiased estimator of a| given by : 

4>lix)dN, = 



"A ^2 



1 /• 1 

^ (/>i(x)d7V, = 5^ 0a(XO'/'a(^p). (3.1) 



Our first approach will consist in constructing estimators of d^(s,5o) = YlxeA '^A based on a 
combination of the Tx's, and in rejecting the null hypothesis when one of these estimators is 
too large. This was already the spirit of Watson's procedure (see |2Z|)- Our second approach is 
related to the test considered in Baraud et al. to detect local alternatives. It will consist in 
considering a set of Ta's and rejecting the null hypothesis directly when one of the Tx's is too 
large. Let us now precisely define both procedures. 



3.1 A first procedure based on model selection 

Assuming that s G L^([0, 1]), a natural idea is to construct a testing procedure from an estimation 
of the squared L^— distance d'^{s,So). 

In order to estimate this functional of s, following the ideas of Laurent [20] and Promont and 
Laurent [13], we introduce embedded finite dimensional linear subspaces of L^([0, 1]). We choose 
here to consider for J > 1 the subspaces Sj generated by the subsets {i;^0)'?^A,A G Aj} of the 
Haar basis defined by fZlJ, with Aj = {(j, k), j G {0, . . . , J - 1}, /c G {0, . . . , 2^ - 1}}. Each 
subspace Sj is called a model. We denote by Dj = 2"^ the dimension of Sj, and by sj the 
orthogonal projection of s onto the model Sj. 

Focusing on one model 5j, we estimate d^(s,5o) = — Oq by the unbiased estimator of 
pjf - «o = EagAj «a given by 

T'j=Yl (3-2) 

AeAj 
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with Tx defined by (|3.ip . The estimator Tj obviously depends on the choice of the model Sj. 
Since we do not want to choose a priori such a model, we consider a collection of models {S'j, J G 
J} where J" is a finite subset of N*, and the corresponding collection of estimators {Tj, J G J}. 
The procedure that we introduce here then consists in rejecting {Hq) "s £ Sq" when there exists 
J in J such that the estimator Tj given by (13. 2|) is too large. 
At this point there are several ways to decide when Tj is too large. 

In all cases, we use the well-known argument that, conditionally on the event "the number of 
points Nl falling into [0, 1] is n", the points of the process obeys the same law as a n-sample 
{Xi, . . . , Xn) with common density s/ ^ s{x)dx. It follows that for all n G N, 

P, {T'j > q'\NL =n)=flj^Yl E MXi)MXi') > q 

y AeAj ii^v=\ 

Under the null hypothesis, the intensity s is constant on [0, 1], and the X/'s are i.i.d., with uniform 
distribution on [0, 1]. This distribution is free from the parameter s. As a consequence, for every 
u g]0, 1[, we can introduce and estimate by Monte Carlo experiments the (1 — u) quantile of the 
distribution of T'j\Nl = n under the null hypothesis, that we denote by qj^\u). 
We now consider the test statistics: 

T« = supfc-,f-)(n'J^-))), (3.3) 

with u^^^ to be correctly chosen. 

Finally, we define the corresponding test function: 

^i^^ = lri^)>o- (3.4) 
And our first test consists in rejecting the null hypothesis (-f^o) when <I>o^ = 1. 
Let us see how we can choose u^J^^^ so that our test has a level a. 



An obvious possibility is to set 



— — - for every J m. J and n in N. 

W I 



This choice corresponds to a Bonferroni procedure and <l?a ^ actually defines a level a. test. Indeed, 
for s G 5o, 



supfc-,r)(n£^)))>0) 



iVi = n P,(7Vl =n) 



nGN JGJ ' ' 

< a. 
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Our choice for u'j^^, inspired by Promont and Laurent [l3], leads to a less conservative procedure. 
It consists in setting 



Hn) _ Wj 
"'J,a — ^ 



sup J u e]0, 1[, sup ( sup (t'j - q'j''\ue-^'')) > 
[ seSo \JeJ ^ ' 



AT^ = n < a 



(3.5) 



where {H^j, J G ^7} is a collection of positive weights such that 

^ e-^^ < 1. 



For the same reason, we still obtain a level a test and by definition, u'j'^-' > ae "^-^ for every n 
in N. 



3.2 A second procedure based on a thresholding approach 

Let us recall here that the squared L^— distance d^(s,5o) between s and the set of constant 
functions Sq is equal to a| and that T\ defined by (|3.ip is an unbiased estimator of a|. 

Based on general thresholding ideas, our second procedure consists in fixing some J > 1 and 
rejecting the null hypothesis (i^o) when there exists A in Aj such that T\ is too large. 
Let us now see what we mean by "Tx is too large". We can still use the fact that 



(Ta >q\NL = n) 



-1 0A(^z)<AA(Xr)>(?) 



and that under the null hypothesis, the X;'s are i.i.d., with uniform distribution on [0,1]. We 
therefore introduce and estimate by Monte Carlo experiments the (1 — u) quantile of the dis- 
tribution of T-x\Nl = n under the null hypothesis, that we denote by q^x ^{u). Notice that for 
A = {j, k) G Aj, q^^\u) does not depend on k. 
We set 

Ti^)=-p(r,-,f^)(ni^-))), (3.6) 

with ti^^^ to be correctly chosen. 
We also define 

^^'^ = lrP)>o- (3.7) 
Our test consists in rejecting the null hypothesis {Hq) when = 1. 

(n) 

Let us now see how we choose u\ ^. An obvious choice corresponding to the Bonferroni procedure 
would be 
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To obtain a less conservative procedure, we prefer setting 

^A^a = — ^ for every A = (j, k) in Aj and n in N, 



23 J 



with 



sup < u g]0, 1[, sup sup r(j- fc) - y, . 1 — ^ 

I se5o V(i,fe)eAj ^ \IJ J 



u 



> 



NL = n\ <a 



When s E 5( 



0, 



(3i 



(3.9) 



sup (t.-,^)(4^)))>0 



sup Tq j,) - 
V {i,fc)GAj V 



(Nl) "a 



> 



> 



NL = n] F{Nl = n) 



< a, 

(2) 

which means that defines a level a test. 
Note that ul"^ > a. 



Comments. 

1. Though the two testing procedures defined by (13.41) and (|3.7p are very different by their 
spirit, they can formally be written in a common way. For any subset A of Aqo, we denote 
by 5a the subspace generated by {(po, cpx, A G A}, by sa the orthogonal projection of s onto 
5a, and we introduce the unbiased estimator = J2xeA'^>^ ~ ckq ~ J2xeA^'x- 

Then our test functions can be written as 

= lr„>o, (3.10) 

where 

Ja — sup Ma ~ '^A,a ) ' 

Aec ^ ^ 
and C is a finite collection of subsets of Aqq. 

Noticing that Tj = T^^, we can easily see that our first test amounts in taking a collec- 
tion C equal to {Aj, J G J'}, and ^Ajo*'^^^ ~ Furthermore, our second 
test amounts in taking a collection C composed of all subsets of Aj, and for A C Aj, 
^'ka^^^ = Z^AeA 9A^^^(^Aa^^)- Indeed, there exists a subset A of Aj such that X^^eA^A > 
EagA <l\'^\'^\a) if and only if there exists A in Aj such that Ta > gi^''^(Mi^'^^). 
Such a common expression will be particularly useful to derive the properties of the tests. 

It also allows us to see our tests as multiple testing procedures. Indeed, we can consider that 
for each A in C, we construct a test rejecting the null hypothesis when T'p^ — t'J^J'^'^'^ > 0. 
We thus obtain a collection of tests and we finally decide to reject the null hypothesis when 
it is rejected for at least one of the tests of the collection. 
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2. Both procedures have a specific interest to prove the optimahty of the lower bounds ob- 
tained in Theorem [D We will actually prove in the next section that the first one achieves 
the lower bounds obtained in case (u) of Theorem [1] (up to a possible logarithmic fac- 
tor) whereas the second one achieves the lower bounds obtained in case (i) of Theorem [H 
However, if we want a procedure that achieves the lower bounds of cases (i) and (ii) simul- 
taneously, we will have to consider the test which consists in mixing the two procedures. 
In this case, we reject the null hypothesis (Hq) when sup{$^\^2> = 1- 



4 Uniform separation rates 

In this section, we evaluate the performances of our new testing procedures from a theoretical 
point of view. More precisely, we prove that our procedures are optimal in the sense that their 
uniform separation rates over Besov bodies are of the same order as the lower bounds for p 
obtained in Section [2l These results justify the construction of our procedures as well as they 
provide the upper bounds needed for the exact evaluation of the minimax separation rates over 
weak and classical Besov bodies in the Poisson framework. 

In the following, the expression C(a, /3, i2, i?', i?", 5, 7, ...) or Ck{a, (3, R, R' , R" , 6,^...) is used to 
denote some constant which only depends on the parameters a, (3, R, R' , R" , 6,^, and which 
may vary from line to line. 

4.1 Uniform separation rates of the first procedure 
4.1.1 The error of second kind 

The aim of the following theorem is to give a condition on the alternative so that our first level 
a test has a prescribed error of second kind. 

Theorem 2. Assume that s S L°°([0, 1]), and that L > 1. Fix some levels a and (3 in ]0, 1[, and 
let be the test function defined by (3^. There exist some positive constants Ci(||s||oo, 
C2{P), C^{a,P), and C4(a) such that when s satisfies 



2 , ^ .„„„ o^vDj , ^ ,^,Dj 



d'{s,So) > inUls - sjW + Ci{\\s\U + C2{P) ^2 



y^^^^^.(x)dx^^^ + — J+C4(a)^^| (4.1 



+ C3(a,/3) 

then 



Comment. Considering here a multiple testing procedure instead of a simple one allows to obtain 
in the right hand side of the inequality (14.11) an infimum over all J in J" at the only price of 
introducing some terms in Wj. These last terms will appear in the following uniform separation 
rates over classical Besov bodies as a In In L factor, which is now known to be the price to pay for 
adaptivity in some classical statistical models. As a consequence, our multiple testing procedure 
is proved to be adaptive in Proposition [Dover classical Besov bodies, which would not occur with 
a simple testing procedure. 
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4.1.2 Uniform separation rates over Besov bodies 

In this section, we evaluate the uniform separation rates p{^a\B2^{R) n L°°(i?"), /?) where p 
is defined by p.ip . and B2^^{R) is any Besov body defined by (12. 2p . 

Let us first notice that the functions oiB2ao{R-) ^re well approximated by their projections onto 
subspaces of the collection {Sj, J G J} considered in our first procedure, in the sense that if 
s G Bl^{R), then 

\\s-sjf < c{5)R^Df^. 

As a consequence we can use Theorem [2] to obtain upper bounds for the uniform separation rates 
of our test. 

We denote by [x\ the integer part of x. 

Proposition 1. Assume that InlnL > 1. Given some levels a and [5 in ]0, 1[, let $0"* defined 
by with J = {!,..., Llog2(LV(ln In L)3)J } and Wj = ln\J\ for every J in J. 

For every 5 > 0, there exists some positive constant C{a, (3, R" ,5) such that when s belongs to 
Bl^{R) n h'^{R") and satisfies 

„^ / VhU^TL \ „2 f (In ^^Lf\ In In L \ 

I 1 J ^ V ) ' 

then 

Fs =o)</3. 

In particular, there exist some positive constants Lo{6) and C{o:, (3, R, R" , 6) such that if L > 
Lo{S), then 

25 

p(^i'\^2,oo(^)nL-(ii"),/3) < C{a,P,R,R!\6) 

Comments. 

1. Our first testing procedure is therefore adaptive: indeed, for large L, it achieves the lower 
bounds for the minimax separation rates over all the spaces B2 ^{R)f^W^{R')r\U^ {R") with 
5 > max(7/2,7/(l + 27)) simultaneously up to a possible InlnL factor (see Theorem [T]). 

However it does not achieve the optimal separation rates obtained in the case where 7 > 
max(2(5, 1/2). In this range of parameters, the regularity in 7 is higher than the regularity 

in 5, meaning that the weak Besov body governs the separation rate. That is the reason 
why we introduced the thresholding type procedure. 

2. The upper bound for the uniform separation rate obtained here is exactly of the same order 
as the (asymptotic) adaptive minimax rate of testing obtained by Ingster [17] in the density 
model, replacing the parameter L of the Poisson model by the number n of observations in 
the density model. In particular the In In L factor is proved to be necessary in the density 
model for adaptive procedures. 



VhilnL 
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3. It is easy to see that Bp f^{R) C ^^(-R) when p > 2. So this result directly leads to upper 

bounds for the uniform separation rates p{^a\Bp,^{R) n h°^{R"), f3) of the test over 
the Besov bodies Bp^^{R) when p > 2. These rates, obtained in the Poisson framework, 
correspond to the ones in some Gaussian models (see Spokoiny |25] for instance) or in the 
density model (see Ingster |17|). 

4. Note that one could also consider some tests based on the Fourier basis as well as the Haar 
basis, as Fromont and Laurent [13] did in the density model. The theoretical results would 
remain unchanged, and the practical performances of the procedure would be better when 
considering smooth alternatives (see Fromont and Laurent |l3j for more details and Section 
E]). We have only considered here tests based on the Haar basis for the sake of simplicity. 

4.2 Uniform separation rates of the second procedure 
4.2.1 The error of second kind 

From the common expression (|3.10ll of the test function for the two procedures, we obtain here 
a result similar to Theorem [2] for the error of second kind of our second test. 

Theorem 3. Assume that s £ L°°([0, 1]), and that L > 1. Fix some levels a and (3 in ]0,1[, 
and let he the test function defined by \3. 7j) . Recall that for any subset A o/Aqo, S\ and sa 
respectively denote the subspace generated by {4>o,4>x,X G A} and the orthogonal projection of s 
onto S\. Denoting by D\ the dimension of Sa, there exist some positive constants Ci(||s||oo, /3), 
C2(/3), Cs{a,P), and C4(a) such that when s satisfies 

d\s, So) > ^mf _ <^ \\s - SAf + Ci{\\s\U(3) ( ^ + JJJ^ j + ^2(0)^ 



^ ^ f . ^Aln(2^J) ,L>A2^1n2(2^J 1 



then 



($(2) = o) < /3. 



4.2.2 Uniform separation rates over Besov bodies 

(2) 

Proposition 2. Assume that InL > 1. Given some levels a and (3 in ]0, 1[, let be the test 
defined by fSJ) with J = [log2(L/ In L)J . 

For every 6 > and 7 > 0, there exists some positive constant C{a, 13, R" ,6,^) such that if s 
belongs to B^2,ociR) ^ W^{R') n ]L°°{R") and satisfies 

then 
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In particular, when 5 > 7/(1+27), there exist some positive constants Lq{5, 7) and C{a, (3, R, R' , R 
such that if L > Lo(5, 7), then 

pi^'i\Bl^iR)nW^iR')nh^iR"),(3) < C{a,P,R,R',R",6,^) 
Comments. 

1. Our second testing procedure is still adaptive: indeed, for large L, it achieves the lower 
bounds for the minimax separation rates over all the spaces B2ao{R) ^ W^{R') D h°°{R") 
with 7/(1 + 27) < 5 < 7/2 simultaneously (see Theorem [T|) . In this case, we also remark 
that these rates are so large that there is no further price to pay for adaptivity in the sense 
that the upper bound does not involve any extra logarithmic factor. To our knowledge, 
this phenomenon is completely new for nonparametric testing procedures. 

2. Our second procedure achieves the lower bounds for the minimax separation rates over all 
the spaces B^^iR) D W^{R') n L°°(i?") with 7/(1 + 27) < 5 < 7/2 simuhaneously, but 
it does not when S > max(7/2,7/(l + 27)). To obtain a test that achieves the minimax 

separation rates in both cases, our two procedures need to be combined. 
4.3 Uniform separation rates of the combined procedure 

Corollary 1. Assume that InlnL > 1. Fix some level a and (3 in ]0, 1[. Let '&^^/2 ^^^^^ 
a/2 test defined by with J = {1,..., [log2{L'^ /{Inln if )\} and Wj = ln\J\ for every J 

in J . Lei ^ Jj^ be the level a/2 test defined by (3. 7\) with J = [log2(-L/lnL)J . We consider 

cI>i=')=SUp{<I>«,cI>i%}. 

(i) For all 6 > and 7 > 0, there exist some positive constants Lq(6) and C{a, /3, R, R" , 6) such 
that if L > Lq{5), then 

24 

Pi'^^a^ , BiooiR) n W^{R') n L-(i?"), /3) < C{a, 13, R, R", 6) 

(a) For all (5,7) such that 7/(27 + I) < 6, there exist some positive constants Lo((5, 7) and 
C{a, (3, R, R' , R" ,6,^) such that if L > Lo((5, 7), then 

p{^^^\BlM n W^{R') n L°°(i?"), /3) < C{a, (3, R, R' , R" , 5, 7) 

Comment. Since 

= 0) < inf {P, (cdW = 0) , P, (c^i^ = 0) } , 

the proof of this result directly comes from Proposition [T] and Proposition [2l 
This final procedure actually matches the lower bounds of Theorem [T] and is consequently adap- 
tive for the whole set of parameters (5, 7) such that 7/(27 + 1) < (5 (up to a InlnL factor 
when 5 > 7/2). This also proves that the lower bounds of Theorem [T] are sharp for this set of 
parameters. 



InL 



1+27 
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5 Simulation study 



We aim in this section at studying the performances of our tests from a practical point of view. 
We consider several intensities s defined on [0,1] such that Jq s{x)dx = 1. N denotes here a 
Poisson process with intensity Ls on [0, 1] with respect to the Lebesgue measure, and the 
distribution of this process. We denote by sq the intensity which is constant (equal to 1) on 
[0, 1]. We choose L = 100 and a level of test a = 0.05. 

Let us now recall that our first procedure may be based on the test statistics 



7c -7 V / 



where qj (u) denotes the (1 — u) quantile of Tj\Nl = n under the hypothesis that s = sq, and 



u'a""^ is chosen such that 



4") =snv\u G]0, 1[,P,„ (sup (t'j - qf\u/\J\)) > 



NL = n\ <a 



The null hypothesis (Hq) "s = sq" is rejected when Ta^^ > 0. 

We choose ^ = {1, . . . , 6}. For 40 < n < 160, we estimate the quantities u'a"'^ and the quan- 
tiles qj"'\u'a"'^ /\J'\) for all J in J". These estimations are based on the simulation of 200000 
independent samples with size n, uniformly distributed on [0,1]. Half of the samples is used 
to estimate the quantiles qj"'\u/\J'\) for u varying on a grid over [0,1], and the other samples 

are used to estimate the probabilities occurring in the definition of Ua"'K Finally, u'a"'^ is esti- 
mated by the largest value on the grid such that these estimated probabilities are smaller than a. 

Let us also recall that our second procedure is based on the test statistics 

AeAj ^ ^ 

where q^^\u) denotes the (1 — n) quantile of Tx\Nl = n under the hypothesis that s = sq. For 
A = (j, k) G Aj, 4"i = u~a^/ {2^ J) with mL"^ defined by dSj]). The null hypothesis {Hq) "s = sq" 

(2) 

is rejected when Ta > 0. 

We choose J = 6. For 40 < n < 160, we estimate the quantities Mq"'* and the quantiles qx^\u^xa) 
for all A G Aj. These estimations are based on the simulation of 200000 independent samples 
with size n, uniformly distributed on [0,1]. Half of the samples is used to estimate the quantiles 
^(Jfc) (^/(^'''^)) ^'^^ ^ varying on a grid over [0, 1], and the other samples are used to estimate 

the probabilities that occur in (13. 9p . Finally, we estimate wi"^ in the same way as in the first 
procedure. 

At this stage, we can estimate the powers of the two tests under various alternatives. The chosen 
alternatives are intensities that have already been studied among others by Reynaud-Bouret and 
Rivoirard [22| . in the estimation problem. Since we are particularly interested in detecting the 
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homogeneity of a Poisson process when the alternatives may be very irregular, we focus on the 
functions defined by: 



si{x) 

S2{x) 



(1 + e)l[0,o.l25[(a;) + (1 - e)l[0.125,0.25[(a;) + 1[0.25,1](^)' 



S3{x) = (1 - e)l[o,i](x) +elY,9j{l + 



■Pj\ 



1[0,1](^) 

0.284 ' 



where 



V = 


[ 0.1 


0.13 


0.15 


0.23 


0.25 


0.4 


0.44 


0.65 


0.76 


0.78 


0.81 


h = 


[ 4 


-4 


3 


-3 


5 


-5 


2 


4 


-4 


2 


-3 


9 = 


[ 4 


5 


3 


4 


5 


4.2 


2.1 


4.3 


3.1 


5.1 


4.2 


w = 


[ 0.005 


0.005 


0.006 


0.01 


0.01 


0.03 


0.01 


0.01 


0.005 


0.008 


0.005 



0<e<l, 0<r/<2, and C2(r/) is such that S2{x)dx = 1. 

These alternatives, for particular values of the parameters, are represented in Figure 1. 
In Figure 2, we represent the histograms of one simulated sample for some of these alternatives 
and for a constant intensity on [0,1]. Note that these histograms are clearly not sufficient to 
separate the alternatives from the null hypothesis. 



s3, ep>silon = 0.5 




u 



Figure 1: Functions si,S2,S3 
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s1 , epsilon=0.7 



s2, eta=1 .5 





Figure 2: Histograms of one simulated Poisson process 

We also consider two monotonous alternatives defined by : 

S4{x) = (1 - e)l[o,o.75[(2;) + (1 +3e)l[o.75,i](2;), 
S5(x) = (1 - e)l[o,i]ix) + s(3x^~H[o,i]{x), 

where < e < 1, and /3 > 1. 

These alternatives, for particular values of the parameters, are represented in Figure 3. 

In Figure 4, we represent the histograms of one simulated sample for some of these alternatives. 



s4, epsilon=0.3 



s5, beta=1 .5, epsilon=0.6 



1 .6 
1 .4 



o.e 
0.4 



1 .6 
1 .4 
1 .2 
1 

0.8 
0.6 
0.4 



Figure 3: Functions S4 and S5 
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Ic3n = 0.3 




lon = 0. 




Figure 4: Histograms of one simulated Poisson process 



For each alternative s, we simulate 20000 Poisson processes with intensity Ls on [0, 1], and we 
estimate the powers of our two tests by : 



1 



20000 



20000 

El 

k=l 



and 



20000 



20000 

El 

k=l 



where Ta^^''^ and T^^'^ are the test statistics 7^^^^ and T^'' computed for the fcth simulated 
Poisson process. 

We compare the obtained estimated powers with the estimated powers of the classical Kol- 
mogorov and Smirnov's test applied to the Poisson process conditionally on the event "the num- 
ber of points of the Poisson process is n". The estimated powers of Kolmogorov and Smirnov's 
test denoted by Pks are also obtained by 20000 simulations of a Poisson process with intensity 
Ls on [0,1]. 

The estimated powers are furthermore compared to the estimated powers of the tests studied in 
practice by the other authors. Such tests are in fact devoted to the particular case of increasing 
alternatives, which may be relevant in reliability contexts involving repairable systems. Bain, 
Engelhardt and Wright [2] and Cohen and Sackrowitz [7] consider in these contexts six well 
known tests. They show that two of these six tests, namely the so-called Laplace and Z tests 
(respectively studied first by Cox |^ and Crow [9]) are preferable to use. 
The Laplace test is based on the statistics 



1=1 



where (Xi, . . . ^Xn^^) are the points of the process, and for every n, q£}^{ci) is the (1 — a) quantile 
of the sum of n independent random variables uniformly distributed on [0, 1]. 
The Z test is based on the statistics 



Nl 



1=1 
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where for every n, (a) is the a quantile of the chi square distribution with 2n degrees of 
freedom. 

Assuming that the intensity s is increasing, the null hypothesis (Hq) "s is constant on [0,1]" is 
rejected when ri^"^ > or T^^^ > 0. 

The readers need to be aware that these tests are especially constructed to detect homogeneity 
against increasing trend, when reading the estimated power tables. 

Let us now present the results we obtained for the different tests. The estimated powers for Pois- 
son processes with intensities Lsi, Ls2, Ls^, Ls^, and Ls^ with various values of the parameters 
are given in the following tables. 



Alternatives si: 








0.5 


0.6 


0.7 


0.8 


0.9 


1 


Pi 


0.05 


0.25 


0.39 


0.56 


0.73 


0.89 


0.98 


P2 


0.05 


0.33 


0.52 


0.72 


0.87 


0.96 


1 


Pks 


0.05 


0.09 


0.13 


0.19 


0.27 


0.37 


0.48 


PLa 


0.05 


0.03 


0.03 


0.03 


0.03 


0.02 


0.02 


Pz 


0.05 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 



Alternatives 



77 





0.5 


1 


1.5 


2 


C2{r,) 


1 


2.27 


3.54 


4.81 


6.08 


Pi 


0.05 


0.61 


0.87 


0.94 


0.97 


7Z 

P2 


0.05 


0.41 


0.64 


0.75 


0.80 


Pks 


0.05 


0.14 


0.25 


0.34 


0.39 


PLa 


0.05 


0.05 


0.06 


0.06 


0.06 


Pz 


0.05 


0.26 


0.39 


0.46 


0.51 



Alternatives S3: 








0.2 


0.3 


0.4 


0.5 


0.6 


Pi 


0.05 


0.28 


0.65 


0.91 


0.99 


1 


P2 


0.05 


0.20 


0.43 


0.71 


0.90 


0.98 


Pks 


0.05 


0.11 


0.21 


0.37 


0.56 


0.76 


PLa 


0.05 


0.01 


0.00 


0.00 


0.00 


0.00 


7! 

Pz 


0.05 


0.02 


0.02 


0.02 


0.02 


0.01 



Alternatives S4: 



e 





0.1 


0.2 


0.3 


0.4 


Pi 


0.05 


0.20 


0.69 


0.97 


1 


P2 


0.05 


0.17 


0.62 


0.95 


1 


Pks 


0.05 


0.26 


0.77 


0.98 


1 


PLa 


0.05 


0.37 


0.82 


0.98 


1 


Pz 


0.05 


0.24 


0.57 


0.85 


0.97 
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Alternatives ^5: 



(/?,£) 


(1.5,0.2) 


(1.5,0.6) 


(1.5,1) 


(2,0.2) 


(2,0.6) 


(2,1) 


Pi 


0.20 


0.49 


0.79 


0.24 


0.62 


1 


P2 


0.18 


0.43 


0.69 


0.24 


0.62 


1 


Pks 


0.22 


0.56 


0.91 


0.24 


0.62 


1 


Pla 


0.24 


0.60 


0.98 


0.24 


0.62 


1 


Pz 


0.24 


0.61 


0.99 


0.24 


0.62 


1 



Comments. 

1. It first emerges from these results that when the alternatives are not increasing, our two 
tests have estimated powers significantly larger than the Laplace and Z tests that are 
designed for increasing alternatives, but also than Kolmogorov and Smirnov's test. Fur- 
thermore, we can not give prior arguments to choose one of our two tests rather than the 
other one in these case. Indeed, we can notice that the first one is more powerful than 
the second one for alternatives S2 which are rather smooth, but also for alternatives S3 
which are very irregular. Thus, in the case of non increasing alternatives such as si, S2 and 
S3, or in the practical situations of our interest such as the study of occurrences on DNA 
sequences where the intensities may have some localized spikes, this should argue in favor 
of the choice of our combined procedure. 

2. As for the increasing alternatives, the specific Laplace and Z tests remain as expected 
the most powerful ones, except for the alternatives S4, that are not as smooth as the S5 
alternatives. Kolmogorov and Smirnov's test is also often more powerful than our tests. 
However, we know that in the case of smooth alternatives, we could probably significantly 
improve the estimated powers of our first test by using the Fourier basis instead of the 
Haar basis. Since our first test is very similar to Fromont and Laurent's p31 one in the 
density model, we refer to this paper for more details. We could also consider a new test 
combining for instance our first test with the Laplace test. 

6 Proofs 

6.1 Proof of Theorem [1] 

Since it is easier to argue in terms of errors of second kind than in terms of minimax separation 
rates directly, we start by defining for all 5 C L^([0, 1]), 

/3(5) = infsupP,($„ = 0), 

where the infimum is taken over all level a tests $q, and by stating a useful and well-known 
lemma. 

Lemma 1. Let r he a positive number, and S , S' he subsets o/L^([0, 1]). 
{i) Iff3{{s e S,d{s,So) > r}) > /?, then 

p{S,a,P) > r. 

{ii) IfS' C S, then f3{S) > /3{S'). 
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The proof of the lemma is straightforward. 

Our aim here is to construct finite sets SM,D,r such that 

SAi,D,r C{s£ Bl^{R)nW.y{R')nh°^{R"),d{s,So) > r}, (6.1) 

and that 

P{SM,D,r) > P, (6.2) 

with r as large as possible. 

These finite sets are based on a family of functions {ifM,i,i S {!,... ,M}} such that for all 
X G [0, 1], ipM,i{x) = (p{Mx — i + 1), where f is a function on [0, 1] such that 

1 

ip{x)dx = 



1 

Lfixfdx = 1 (6.3) 



VxG [0,l],|v;(x)| <p. 
For r > 0, and D < M, we introduce the set 



Sm. 



.D,r = < S?,A,r = /Ol[0,l] + ^ ^i^i^M^i, C S {-1, A G {0, = ^ \ ' 

I i=l i=l ) 

(6.4) 

As a first step, we notice that the functions s^,A,r's are positive as soon as < D /M and that 

for every sj^A.r £ SM,D,r, (^Isg.A.r, 5o)^ = ||s^,A,r - pl[o,i]lP = (see ([Oil ). 

As a second step, we want to find which positive r leads to (3{SM,D,r) > /?• 

Let us recall a fundamental lemma which can be found in Ingster [16] or Baraud p] for other 

frameworks. 

Lemma 2. Let v he a probability measure on SM,D,r and let a ~ u. Let be the distribution 
of a point process N such that the conditional distribution of N given that a = s is a Poisson 
process with intensity s. Let Pq be the distribution of a Poisson process with constant intensity 
given by pi [0,1], and Eq denote the expectation with respect to Pq. Let Ly be the likelihood ratio 
Ly = dP^/dp'o- Then 

P{SM,D,r) > 1 - a - ^ (MLliN)] - 1)'/' . 
Proof. The proof is obtained by rather straightforward computations. One has 

P{SM,D,r) > inf / P.($a = 0)diy{s) 

*Q J 

> 1 - sup / Fs{<^a = l)du{s) 

> 1 - sup [\F,{^a = 1) - Po(^a = 1)1 + |IPo(^a = 1)|] 



> 1 - a 



II TV, 
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where \\-\\tv corresponds to the total variation norm. Hence, 



1, 



> l-a--Eo 



P{SM,D,r 

> l-a-^Eo[\U{N)-l 
But Eo[L,{N)] = 1. So (3{SM,D,r) > 1 - a - {Eo[LI{N)] - 1)^/^/2. 



1/2 



Regarding Lemma[2l we still have to find a distribution v and r such that E,q[Li,{N)] < 1 + 4(1 — 
a — (3)"^ which implies that I3{Sm,d,t) > 

Let ^ = (^1,... ,Cm) be a random vector, such that the ^j's are i.i.d. Rademacher variables, 
taking the values +1 and —1 with probability 1/2. Let A = (Ai, . . . , Ajv/) be a random vector, 
independent of ^ and defined by Aj = Ijgx, where X is a set of D indices drawn at random from 
{1, . . . , M} without replacement. 

Then the random function s^,A,r = P'^[{),i\~'rf'\J^YliLi ^i£,i^M,i belongs to SM,D,r, which allows 
to take its distribution as z^. 

Let us denote by the expectation with respect to the variable ^ and by Ej the expectation 
with respect to the random set I defined above. By definition, Ly = J dFs/dFQdi'{s). Hence 



Ly{N) = EjEg exp ( \n{s^^^^r{x) / p)dK 



This can be rewritten as 



Ly{N) = ExE. 



■ M 

JJexp 



i=l 



i — 1 i 
M ' M 



dN^ 



E7 



Ha 



where 



exp 



In 1 + ri 



i-l J_] 
M ' M J 



D p 



dN, + 



exp 



In 1 — ri 



'MlpmAx) 



dN,, 



i-\ J_] 

M ' M \ 



D P ^ 

Let I' be a random set of indices with the same distribution as I and independent of I. Then, 



Eo[Ll{N)] = EoEjEi: 



iei iei' 

= EoExEr[ n n n ^^'i 

ieJ\i' ieT\i ieinx' 

But under the distribution Pq, the variables Aj's are mutually independent since they only depend 
on the integrals of the Poisson process on intervals with disjoint support. Consequently, 



Eo[Ll{N)]=EjEr 



II Eo[Ai] Yl Eo[A,] Yl Eo[Aj 
iei\i' iei'\i ieini' 



(6.5) 



We now need to compute Eo[^i] and Eo[A?], and we use the following lemma. 
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Lemma 3. Let f be a function on [0, 1]. Then with the above notations, 



En 



exp 



AI 'Ml 



f{x)dN^ 



exp 



i-l J_] 

M ' M J 



(exp(/(x)) - 1) pLdx 



Proof. When has the constant intensity pl[o,i], we know that conditionally on the event 
"the number of points A^M,j = (] if]) falling into ] j^] is n", the points of the pro- 
cess Xi,. . . jX^Mi in ] ^) Ji] obey the same law as a n-sample with uniform distribution on 
] Jl] ■ Then, one can easily see that 



En 



exp 



i-l 

M 'Mi 



f{x)dN, 



En 



En 



En 



exp Yl fi^i 



1=1 



NM,i 



n E[exp(/(XO)|iVA./,, 



1=1 



exp ( A^A/,i In I / exp(/(x))M(ix 

aF' lu] J 



Under Pq, -^A/,i has a Poisson distribution with parameter pL/M ^ therefore, 



En 



exp Nma In 



i-l i 1 
M 'Mi 



ex.p{f{x))Mdx 



pL 



exp I ^ I / ^ exp(/(x))Mdx - 1 



This concludes the proof. 
From Lemma M and (16. 3|) , one has that 

1 



MM = -exp 

2 \/tiij_\Vi' p 

" J M ' M J \ 



1 

+ o exp 

2 \ /I izii J_ 

1. 



-r^^^^^\pLdx 



Moreover, 



Eo[Aj] = -Eo 



4 



exp 



+^Eo 



+^Eo 



exp 



21n 1 + diV. 
'-i ^1 \ ^ D p 



21n|l-rW^WMUAr. 



z -l J_] 
M ' M J 



exp I / In 

M ' M J 



1 — r — — 



D p 



2 , dN^ 
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Using Lemma [3] and (|6.3p again, we finally obtain that 



Eo[Af] = cosh 



- 



Hence, equation (I6.5|l gives 



cosh 



7d) 



exp ( |X n X'l In cosh 



pD 



For fixed X, |XnX'| is an hypergeometric variable with parameters {M,D,D/M). Hence, we 
know from Aldous [l| p. 173, that there exists a binomial variable B with parameter {D,D/M) 
such that Ej/ [i?||XnX'|] = |XnX'|. By Jensen's inequality, we obtain that 



Eo[Ll{N)]<EiEr 



exp ( B In cosh ( 



Setting B = Yld=i Bi where the Sj's are independent random Bernoulli variables with parameter 
D /M , we easily obtain that 



Eo[Li{N)] < exp ( L» In ( 1 + ^ ( cosh 



M 



-1 



From equation (16. 6|) and Lemma [21 we see that if 



exp D In 1 + 



D 

M 



^cosh ^ 



pDj 



1 <l + 4(l-a-/?)^ 



(6.6) 



then (i{SM,D,r) > P- 

Following Baraud's idea \2l and setting c = 1 + 4(1 — a — since the function cosh is increasing 
on [0, +oo[, we have that if 



n pD I M I M (M 

r < ^ In 1 + — r In c + W 2— r In c + I — In c 



L 



Z)2 



Z)2 



L>2 



then 



cosh I — — - 
\pD 



1 / M , / M , / M , 

^ ^ 2 b^'"^+vV^^^+U^"^ 



< 



Mine 
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Hence 



exp(i?ln(l + ^(cosh(^)-l 



< exp L» In 1 + 



< c, 



In c 

IT 



and (5{Sm,d,t) > /3- As a conclusion, we obtain the following result, where the second part of 
the Proposition comes from a direct computation (see Baraud [3] for further details). 

Proposition 3. Let c = 1 + 4(1 — a — (3)'^ and SM,D,r be the finite set defined by {6.4^ . If 



D 



pD 



M 



M 



r < —— and r < —— In 1 + —77 In c + \ \ 1—^ In c + ( —77 In c 



M 



L 



L>2 



M 



L>2 



(6.7) 



then SM,D,r C {s,s > 0,d{s,So) = r} and f3{SM,D,r) > P- 
//« + /? <' 0.59 and 



2 D 

r < ITT /\ 
- M 



pD I M M 



then SM,D,r C {s,s> 0,d{s,So) = r} and /3{Sm,d,t) > P- 

As a third step, we are now in position to find some r (as large as possible) such that Sm,d,t C 
{s £ Bi^{R) n Wj{R') n h'^{R"),d{s, So) > r} and that f5{SM,D,r) > P- 



Let us consider the set SM,D,r defined by (|6.4p with ip = l[o,i/2[ ~ l[i/2.i[i M = 2"^ and p = 1- 



Let s G SM,D,n then s can be rewritten as s = ao'/'o + SjeN X^fcLo^ o:(j,k)4'{j,k)i with 



"0 = 1, ay,k) =Oif j J, afj,^^ = -jj^k+i for k = Q, . . .2'^ - I. 

Since Xlfcli ^k = D, the condition < R^2^'^-^^ ensures that SM,D,r C Bi ,^{R). 
Let us define, for all t > 0, 

2^ 2 

H{t)=Y.''-/^klr^ 



k=l 



D " i^Afc<i 



In order to ensure that s belongs to W-y{R'), the function H has to satisfy 

yt > 0,H{t) < Rl'^t^. 



Note that 



2 / 2 \ 2 

H{t) = for f < ^ and H{t) = H y — \ for t > ^. 



Hence, we only need to have that 
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which is equivalent to 

Moreover, if < D/M, then ||s||oo < 2. Hence when R" > 2, the condition 

r^<^A R^M-^' A i^'2(i+27)^-27 (g_8) 

ensures that SM,D,r C Bl^{R) n Wj{R')n'L°°{R"). Prom Proposition El we can conclude that 
when i?" > 2 and a + /3 < 0.59, if 

< ^ A (^f In (^1 + ^ V ) ) A R^M~^^ A R'^i^^^^^ D-^^ , (6.9) 

then (|6.ip and (|6.2p are both satisfied. 

We now consider several cases, that are represented on the following figure. 




In the following of this proof, C will denote a positive constant that may depend on a, /3, R, R' , R", 6, 7, 
and that may vary from one line to another. 



Case 1. If (5 < 7/2, and 6 > 7/(1 + 27), we set 



D 



L \ 1+27 



InL 



25 



and 



M = 2"', with J ■ 



log2 



j/S 

L \ 1+2^ 



InL 



+ 1. 



We first check that D < M for L large enough since S < j. 
Then, 



and 



Finally, since 



then 



and 



-27 



M 7/'5-2 

> (L/lnL) 1+27 ^L^+oo +00 when 7 > 25, 



D M M \ -27 

+ >C(L/lnL)TT^, 



1-7/^ -27 
— > C(L/lnL) 1+27 > C(L/lnL)i+27 for L large enough. 



Case 2. If 7 > 1/2 and 5 < 7/(1 + 27), one chooses 



D 



L \ 1+27 



InL 



and 

M = 2-^ with J = [log2(L/lnL)J + 1. 
We first check that D < M for L large enough since 7 > 0. Then, 

27 



and 

Since moreover 



R^M-^^ > CiL/liiL)-'^" > C{L/\nL)^. 



25 



-27 



and 



M 1 2_ 

— > (L/lnL)' 1+27 ^L_^+oo +00 when 7 > 1/2, 

D ( M [W\ -27 
+ >C(L/lnL)iT^, 



Z) -27 

— > C(L/ InL) 1+27 for L large enough. 



Case 3. If 5 < 7 < 25, and 5 > 7/(1 + 27), one chooses 

M = 2^ with J = log2(L2/(i+4^)) 



+ 1, 



26 



and 

D = 

With such a choice, one has that D < M and 



^/2(l+27)£,-27>c^T^^ 

Furthermore 

^ ~ ^^^^^ when (5 > 7/2 and 1 when S = 

Hence, 

L> / M fM'\ _ ^ 

when 5 > 7/2, and 

D r/ 1 -4^ 

iW^/^-^ > 

M 

when 5 > 7/(1 + 27). 

Case 4. If 7 < 5, one chooses M = D = 2-^ with J = [log2(L2/(i+45))J + 
With such a choice, 

^/2(1+2t)^-27 > ^/2(l+27)^-25 > CL^^, 
M 



and 

Moreover 



so 



and 



D -is 
— = 1 > Li+«, 
M 

for L large enough. 

Case 5. If 7 < 1/2 and 5 < 7/(1 + 27), one takes 

M = 2'' with J = [Iog2 LJ + 1, 

and 

We first notice that D < M. Then, 

^/2(l+27)^-27>(^^I^^ 



27 



and 



Moreover 



Hence, 



and 



M 27-1 

_ -^L^^^ when 7 < 1/2, and 1 when 7 = 1/2. 



D -27 -I 

— > CL^ > CL-'/' for L large enouen. 

M ~ ~ 6 6 



This concludes the proof of Theorem [H 



6.2 Proofs of Theorem H and Theorem SI 
6.2.1 Preliminary results 

We consider here the general test function $q = lra>o, defined by (|3.1Q|) . where 

la — sup I - tj^ 

Aec ^ 

Tj^ = Y1x£a'^>^^ ^^'^ C is a finite collection of subsets of Aqo- The collection C and the quantile 
^'ao^^ will be chosen to fit our two procedures respectively. 
We begin to prove the following result. 

Theorem 4. Let s G L°°([0, 1]), and fix a and (3 in ]0, 1[. Assume that there exists some positive 
quantity ^A,a,/3 such that 

We recall that Da denotes the dimension of S\ and we set E\ = '^j/y k)eA ^'^ ■ 
There exist some positive constants Ci(||,s||oo, /?) and C2{(3) such that when s satisfies 



d\s,So) > mf <j lis - SAir +Ci(||s||oo,/3) + + C2(/3)|^ + ^A,a,/3 



(6.10) 



then 

Ps (^a = 0) < /?. 

Proof. Let a and /3 in ]0, 1[, and s be a fixed intensity. 



($„ = o) = Ps(r„<o) 

vAGC,n'<t?;^-) 



< infP.(T;:<tr^) 
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For every A in C, we can write Tj^ in the following way : 



AeA 



[0,1] / J%i\ 



L2 



AeA 



0a(2;) ((iA'x - +2/ (j)x{x)dN^ (j)x{x)s{x)Ldx 

,1] / ^[0,1] ^[0,1] 



L2 



AeA 



[0,1] 



0A(x)s(a;)Ld2; ) + / (j)l{x)dN, 
[0,1] 



By setting 



AeA 



,^a(x) {dN, - s{x)Ldx) - / <Pl{x)dN, 
[0,1] / ^[0,1] 



and 



L 



[0,1] 



(sa(x) — aQ(t)Q{x)) {dNx — s{x)Ldx) 



we obtain the following decomposition : 

T'jl = UK + VA + \\sAf -al 
Since d^(s,5o) = ||s — sa|P + ||sAp — follows that 

r;' = Uk + Vk + d\s,So) - \\s - sAf. 

Hence, 

P. (^a = 0) < inf Fs (Ua + Va + d\s, So) < \\s - SAf + t^^^^^ 
AeC V ' 



(6.11) 



The aim of the following lemmas is to define positive quantities A^j^^^ and A^^^^, such that 

P. (^A < < f , 

P. (Va < -A?^.) < i 



A,PJ ^ 3 ■ 



Using (16.111) and assuming that 



we then obtain that as soon as there exists A in C such that 

j2/„ c ^ ^ ll„ „ . I|2^ .(1) , 4(2) 



then 



d'{s, So) > lis - SAir + ^A,/? + ^A /3 + 
Ps (^>a = 0) < /?. 



(6.12) 
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Lemma 4. There exists some positive constant C such that for all K ^ C and for all a; > 0, 



^ -t/A > C \\sUl^^+\\sU^x + + j < 2.77e- 

Proof. Let us first notice that 



Setting 



L2 



AeA 



[0,1] 



'[0,1] 

(t)x{x)dN^ 



ct)x{x)dN^ 



[0,1] 



ct>l{x)dN^ 

'[0,1] 

4>\{x)s{x)Ldx + 



[0,1] 



(f)\{x)s{x)Ldx 



1 r f r 



L2 

+ 1 

2 

l2 



(i)x{x)dN.^ 



[0,1] 



cl)\{x)s{x)Ldx 



[0,1] 



(j)\{x)s{x)Ldx 



E 



aga 



1 ry 



JO 



(/.A(x)</.A(y)c?7Va;C^iVy 



1 



+ 



JO 

1 /-y- 



1 /-i 



<l)x{x)(l)x{y)dN^s{y)Ldy + / / (t)x{x)(l)\{y)dN^,s{y)Ldy 



Jy 



Jo 



4>x{x)4>x{y)s{x)Ldxs{y)Ldy 
1 /■!/- 



aga L 
^1 r?/ 



JO 



JO 

1 ry~ 



1 fX 



JO 



JO 

-y 

Aga 



0A(2;)0A(y)s(a;)-^^'^a;s(y)L(iy 
1 ry~ 

<Px{x)(t>x{y){dN^ - s{x)Ldx){dNy - s{y)Ldy) 



Jo 



2 

HA{x,y) = j2'^'Pxix)4>x{y), 



(6.13) 



aga 



we deduce from Theorem 4.2 in Houdre and Reynaud-Bouret [15| that there exists some absolute 
constant k > such that for all x > 0, 

Ua> h(Ai^/^ + A2X + ^3x^/2 + A^x'^) ) < 2.11 e-"" , 
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where 



Ha{x, y)s{x)Ldxs{y)Ldy, 

sup / a{x) i b{y)HA{x,y)s{y)Ldy ] s{x)Ldx, 

x)Ldx=J b'^{x)s{x)Ldx=lJO \Jx 

sup / H\{x,y)s{x)Ldx, 
GfO,ll -'0 



Ao = 



^3 



ye [0,1]' 
sup \HA{x,y)\. 

x,y£[0,l] 



Let us now evaluate ^1, A2, and A\ for every A G C 
To give an upper bound for A\, we notice that 



Al<\\s\l^L^ [ [ Hl{x,y)dxdy. 

Jo Jo 

Since {(f)\, A G A} is an orthonormal basis on [0, 1], one has 

■ /' I'yZM^fMyfdxdy 

Jo Jo 



Al < ^2 



< 



M\s\\lD 



A 



L2 



For A2, we use Cauchy-Schwarz inequahty to see that 

"1 / rl 



A2 < sup 

bJb'^ix)s{x)Ldx=l 







J Ky)H\{x,y)s{y)Ldy^ s{x)Ldx 



1/2 



and 



j\{y)HA{x,y)s{y)Ldy^ < (^j\''{y)s\y)Ldy^ (^J^ Hl{x,y)Ldy 



This imphes that 



Ao < L 



J \\4oo\^J Hl{x,y)dyj s{x)dx 

nH\{x,y)dxdy 



1/2 



< m 

Since {^a, A G A} is an orthonormal basis on [0, 1], one has 



A2 < -^||s||oo 



r^4>l{x)cl^liy)dxdy 

•^0 -^0 AeA 



1/2 



< |lHloo<. 
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As for ^3, we can prove in the same way that 

4||s 



3 ^ —FY- ™p ^^liy)- 



Moreover, for any y fixed in [0, 1], 



This imphes that 
Furthermore, for x,y in [0,1], 



(i,fc)eA j/U,k)eA 
< Ea. 

.2 ^ 4||s||oo^A 



|i?A(x,y)| 



2 

l2 



(i,fc)eA 



i/(i,fe)eA 



< 



2Ea 
L2 ■ 



Finally, A4 < 2E\/L?' , and this concludes the proof of Lemma [H 

By taking x = ln(8.31//3) in Lemma HJ we obtain that a possible value for A^^^ is 



A^"-^ - C WsW - 



-|^21n(8.31//3) + ^^^(ln(8.31//3))3/2 + |^(ln(8.31//?))^ 



where C is an absolute positive constant. We now use the following lemma, which derives from 
an analogue of Bennett's inequality (see proposition 7 of Reynaud-Bouret [2T], for instance). 

Lemma 5. There exists some positive constant C such that for all x > 0, 



■ s I -Va > - ao(, 



f-l\\s-SAf + ^^x]<e'^. 



Proof. Recall that 



L 



[0,1] 



(sa(x) - ao(f)o{x)) {dNx - s{x)Ldx) . 



Using proposition 7 of Reynaud-Bouret pi], we easily obtain that for all x > 0, 



-Vk>2\2x 



L 



II J. 112 , ^i'^A - ao<} 
|SA - ao'Poi H 



-X < e 
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First note that 

||SA - ao4>o\\oo < \\s\\oo- 
By using the elementary inequahty 2ab < + 26^, we obtain that 



2\/ 2x^-^^^^\\sA - ao4>of < ;^||sA - ao0o||^ + 4x- 



L " ^ " - 2" ^ " L 

1 11 / ii2 II II 2 A II ^ II cxD 

We deduce that there exists C > such that for all a; > 0, 

-Va > ^||. - aoM' -lis- sj,f + 9M21^ < e-. 



(2) 

By taking x = ln(3/^) in Lemma [5l we obtain that a possible value for is 

= - ^oM' -lis- SAf + ^ ln(3//3). 

Replacing ^j^^^^ and /? ™ (|6.12p by the possible values obtained above finally leads to the 
result of Theorem [H ■ 



We now prove the following lemma that will provide an upper bound for the quantity ^A,a,/3 
occurring in Theorem [H 

Lemma 6. Let Xi, . . . ,Xn be i.i.d. uniformly distributed on [0, 1]. For n G N and A C Aqo, let 

n 

Tin = E E UXMxiX,). 
AeA/7^/'=i 

Let Da denote the dimension of Sa and Ea = fc)gA ^-^ • There exists some absolute constant 

C > such that for all 2; > 0, 

P {T'ir. > ^ + X + f^)) < 2.77e-. (6.14) 

Proof. If n G {0, 1} , T'^ ^ = Q hence (I6.14|) holds. Since for all A G Aqo, 4>\ is orthonormal to 
4)0 = l[o,i], it follows that the variables (j)\{Xi) are centered and we can apply Theorem 3.4 in 
Houdre and Reynaud-Bouret ^5]. We now set HA{x,y) = Ylx^A 'l^>^(^)^>^(y) / ^'^ ■ obtain that 
there exists some absolute constant C > such that for all x > 0, 

P (rl^ > C [Ai^ + A2X + isx^/^ + Mx^)) < 2.77e-^ 

where 
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H\{Xi,X2) 

n l-l 

^^/7A(Xi,X2)a;(Xi)A'(^2) 



E 



A2 = sup 

j=i /'=i 

Al = n sup / Hl{x,y)dx, 

j/G[0,l] Jo 

= sup \HAix,y)\. 

x,yg[0,l] 



,E 



< 1,E 



ii=i 



< 1 



To evaluate Ai, A2, A3, A4, we use arguments similar to the ones used in the proof of Lemma IH 
Since {(px, A G A} is an orthonormal basis on [0, 1], 



2 /•! /■! 



Jo 



AeA 



< 



n^DA 



Let («!,..., On) and . . . , such that E Yl?=i'^'ii-^i) 
Then 



< 1 and E 



< 1. 



" n l-l 
J=l i'=l 

= EE/ / HA{x,y)ai{x)(3i'{y)dxdy 



1=1 i'=i 



JO 



n 



TsXiyiy] / 0A(2;)ai(2;)d2; / (i)x{y)[ii>{y)dy. 



By using Cauchy-Schwarz inequality, we obtain 

n l~l 



E 



< 



/7a(Xi , X2)a/ (Xi) A' (^2) 

_l=l V=i 

^ n l-l V y „i 

^^Yl Y{ <P\{x)ai{x)dx 
1=1 v=i Laga ^ ° 



1/2 



y^( f\x{y)Pi'{y)dy 



1/2 



One has for all g G L^([0, 1]), X]AeA(/ '/'A^)^ < / 5^- As a consequence, 



1/2 n 



/'=i 



1/2 



< 



n 
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We evaluate and in the same way as and A4 in the proof of Lemma [H We obtain that 

^3 < ^4 ' 

and 

Ai < E^/L\ 

Finally, we proved that there exists some absolute constant C > such that 



Since 



'n n 

we can simplify the above inequality : there exists some constant C > such that 



nV^A 3/2 , , 2 
2 — —X ' < X -\ X , 



> {yo^ + X + ^x^^ ^ < 2.77e- 



for all X > 0. This concludes the proof of Lemma [H ■ 
We are now in position to prove Theorem [2] and Theorem [3l 

6.2.2 Proof of Theorem M 

Recall that the test function defined by (|3.4p is of the same form as the test function (IS.lOp of 



Theorem m with C = {Aj, J G J'}, and = Q'f^' {u'j^J), where q^^' {u) denotes the (1 

quantile of T'p^^ ^. Since Uj]l defined by (13. 5|) satisfies u'f'^ > ae~'^-' for all n, one has that for 
all n, 

In order to use Theorem IH we then need to find some positive quantity ^j,o,/3 such that 

P,(,f-)(ae-^0>^w)<f- (6.15) 

Let us first give an upper bound for qj"'\ae~^-') for all n in N. We apply (I6.14p with A = Aj 
(note that L>Aj = -Eaj = Dj) and with x = ln(2.77/a) + Wj. There exists some absolute 
constant C > such that 

Ve-'^O < C-^ {^Dj {\n{2.7V a) + Wj) + ln(2.77/a) +Wj + -^(ln(2.77/a) + Wjf^ . 
This allows us to obtain some Aj^a,p such that (|6.15p holds. It actually gives that 
gJ^^)(ae-^0 < {^Dj (ln(2.77/a) + Wj) + ln(2.77/a) + Wj] + ^{\n{2.n/a) + Wjf 



Now, from Bernstein's inequality, we deduce that for all n > 0, 



\Nl> I s{x)Ldx + . 2[ s{x)Ldxu + - 



3U <.-". 
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Hence a possible value for Axj,a,i3 is 



frr, -,-1 s(x)Ldx + ln( 3/ B) , , x n. 

12 (V^J (ln(2.77/a) + Wj) + ln(2.77/a) + T^^j + -^(ln(2.77/a)+Tyj)^ 

Using Theorem H] finally leads to the result of Theorem [2l 



6.2.3 Proof of Theorem [3] 

Recall here that the test function defined by (13. 7|) is of the same form as the test function (I3.1QP 
of Theoremllwith C = {A, A C Aj}, and = ^AeA (u^S^ l{VJ)^ , where q-^^ {u) denotes 
the (1 — u) quantile of Tx conditionally on the event = n under the null hypothesis (-ffo) and 
u'^^ defined by (13. 9|) satisfies ui"^ > a for all n. 



Hence, we can prove Theorem [3] by using Theorem [4] and some positive quantity ^A,a,/3 such that 



Following the same lines of proof as in the previous section, let us first give an upper bound for 
«S) («/(2^'^"))- 

Notice that (t^^^^jyU) is the (1 — u) quantile of the variable Tp ^ with F = {(j, A;)}. Since Dy = 1 

and Er = 2^ , the inequality (I6.14p implies that there exists some constant C > such that for 
all X > 0, 

F(r^„>C^(v; + . + 2^^))<277e-. 
Taking x = ln(2.77) + ln{2^ J/a) in this inequality leads to the conclusion that : 

q'Sl) (^) < (^(ln(2.77) +ln(2^ j7a)) +ln(2.77) +ln(2V7a))+|i(ln(2.77)+ln(2V7a))^ 
From Bernstein's inequality, we deduce that a possible value for A\^a,p is 

C \^ J2 f7(ln(2.77)+ln(2iJ/a))+ln(2.77)+ln(2V/a)j 



(j,fc)eA 



+ |^(ln(2.77) +ln(2V/a))2 I, (6.16) 



for some positive constant C. 

Since |A| = Da — 1 and Ea < 2"^, we obtain the result of Theorem [3l 
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6.3 Proof of Proposition [T] 

Let us assume that s belongs to S2oo(^) ^ L.°°{R"). We need to find an upper bound for the 
quantity 



ml {\\s-sj\\' + Ci{\\s\\oo,^)^ + C2iP)-^ + C3ia, 13) / s{x)dx 



L2 



[0,1] 



ypjwj wj 

L L 



+ CJa) 



DjW] 



L2 



in Theorem [2j 

We have already noticed that when s belongs to B2^^{R), for all J > 1, 

Moreover, the constant Ci(||s||oo, /3) can be replaced by Ci{R",(3), so we only need to find an 
upper bound for 



2^-25 , , Dj , VDJW] , Wj , DjW] 



C{a, P, R", 6) inf { R^Df^ + + JJ + 



L 



+ -+ L2 



Taking Wj = ln\J\ = In [log2(LV(lnlnL)3)J , with InlnL > 1 leads to Wj < 2.06 In In L, so 



Jej 



C{a, (3, R", 6) inf { R^Df^ + + + 



Dj Dj , VDJW] , Wj , DjWj 



L L2 
<C'(a,P,R",6) ( inf \r^D-'^^ + 
Since for all J in J, Dj < L^/{lnlnLf, 



L 



+ -+ L2 



^/Dj InlnL ^ L>j(lnlnL)2 



L2 



+ 



In InL 



C'{a, P, R",5) jn^ li^^i^y^-^ + 



2„_25 V-Dj InlnL L>j(lnlnL)2| InlnL 



L 



+ 



L2 



+ 



L 



< /3, i?", S) ( int ^ i?^I)7^^ + ^5Z^ I ^ In In L 



L 



L 



We have that i222;j2<5 ^ In In L/L if and only if J > logs ((i^^LVlnlnL) 
we introduce 

and we distinguish three cases. 

When 1 < J* < Llog2(LV(lnlnL)3)J , then J* belongs to J and 



1/(1+45) 



Hence, 



inf lR'Df'+^^^^^™^ 



L 



< ^2^j^2. + y^i, InlnL 



< (1 + V2)R^ 



L 

Vln In L 
L 



4,5 
4(5+1 
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When J* > Llog2(LV(lnlnL)3)J , this means that for all J in J, y/Dj In In L/L < R^Dj^^. By 
taking J* = |_log2(L^/(ln In L)^)J , we obtain that 

in^ i^R'Df' + ^^Z^l < 2R'DJ^' < 2^'^^R' {^^^^^ 

Finally, when J* < 1, then for all J in J , R^DJ^^ < \/Dj In In L/L, so by taking J* = 1, we 
obtain that 



, f o 9^ V^J nlnLl V21n nL 
Jej 1 J ^ 



This ends the proof. 

6.4 Proof of Proposition [2] 

Let us assume that s belongs to n n L°°(i?")- We now need to find an adequate 

upper bound for 

^mf^_ jlk - ^aII' + C,{\\s\U(3) + f^) + C2{f3)^ 

+ C3{a,f3)J s{x)dx hC4(a) j-^ > 

in Theorem [3l 

As in the proof of Proposition [H the constant Ci(||s||oo, /5) can be replaced by Ci{R",f3). More- 
over, with the choice J = [log2(L/ In L)J , we have that 

2J/2 I 
< 



L3/2 - LVh^L' 
2^ 1 



< 



L2 - LlnL' 

and In (2'^ J) < InL. So we only need to find an upper bound for 



L 



Let us introduce for all integer D <2^ the subset A/j of Aj such that the elements of {a\, A G 
A/)} are the {D — 1) largest elements in {oa, A G Aj}. 
We can notice that 



II ||2 _ II _||2 , II _ ||2 



On the one hand, since s belongs to i32oo(-R), 

-25 



\s - sjf < C{5)R^ 



L 
h^L 



38 



On the other hand, since s belongs to W^{R'), then for all t > 0, 



Yl Yl ^i«o-.<c)i>i - 



jeN fe=o 



< 



2^-1 

YYYl l|2'<|ao,,)|<f2'+i 
jGN k=0 leN 

ieN jeN fe=o ^ 2 



l"(i,fc)l 



I«0,fe)l<f2'+i 



-2/ 2^-1 



ZeN jgN fc=0 

12 



27 

2i \ 1+27 



Taking t such that C{'y)R t 1+27 = Z) in the above inequality proves that all the coefficients of 
sj — are smaller than t/2 and 



J-l 2J-1 



j=0 k=0 
< C(7)ii:'2+47£)-27. 



Hence, 



<C(a,/3,i?",5,7) ( inf <| i?'2+47£,-27 + 

\ 1<D<2J I ^ 



L»lnL 



-25 



+ 



InLy lVI^ 



We have that i?'2+47£)-27 < DlnL/L if and only if L» > i?'^ (L/lnL)^/(^+^^\ Hence, we 
introduce 

1 

L \ 1+27 



72 



InL 



+ 1, 



and we distinguish two cases. 

When 1 < < 2^^, we clearly obtain that 



inf \ i?'2+47^-27 + I < i?'2+47^-27 + £l^_ 



On the one hand, when D^, >2, this leads to 
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1<D<2J L I L 



L \ 1+27 



On the other hand, when = 1, since one has 

inf Ir''-^'^D-''' + ^]<2^^. 
1<D<2J [ L j L 

Now, let us consider the case where > 2'^. This means that for all D such that 1 < -D < 2"^, 
DlnL/L < /2'2+47£)-27. By taking D* = 2\ we obtain that 



1<D<2 



T \ -27 
21nL 



This concludes the proof of Proposition [2l 
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